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O ' ABSTRACT 

u: 

Context. This paper presents a reconstruction of the gravitational lens SL2S02176-0513 
using the singular perturbative method presented in Alard 2007, MNRAS Letters, 382, 
58 and Alard, C, 2008, MNRAS, 388, 375. 

Aims. The ability of the perturbative method to separate the inner and outer contri- 

butions of the potential in gravitational lenses is tested using SL2S02176-0513. In this 

lens, the gravitational field of the central galaxy is dominated by a nearby group of 

' ^ , galaxies located at a distance of a few critical radius. 

Methods. The perturbative functionals are re-constructed using local polynomials. The 

polynomial interpolation is smoothed using Fourier series, and numerically fitted to 

HST data using a non-linear minimization procedure. The potential inside and outside 

the critical circle is derived from the reconstruction of the perturbative fields. 

Results. The inner and outer potential contours are very different. The inner contours 

are consistent with the central galaxy, while the outer contours are fully consistent 

^ . with the perturbation introduced by the group of galaxies. 

•rH . 

Conclusions. The ability of the perturbative method to separate the inner and outer 
contribution is confirmed, and indicates that in the perturbative approach the field 
of the central deflector can be separated from outer perturbations. The separation of 
the inner and outer contribution is especially important for the study of the shape of 
dark matter halo's as well as for the statistical analysis of the effect of dark matter 
substructures. 

Key words, gravitational lensing - cosmology:dark matter 



1. Introduction. 

Gravitational lenses are a precious astrophysical tools to probe the mass distribution of dark 
halo's. The recent results obtained by Clowe etal (2006) are a good example of the ability 
of gravitational lensing to probe complex mass distributions. In this context the recently 
perturbative approach (Alard 2007), and its application to the reconstruction of strong 
gravitational lenses (Alard 2008, 2009), is of particular interest since this approach relates 
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directly the properties of the images to the lens potential. The lens SL2S02176-0513 is the 
second lens to be reconstructed using the perturbative approach. The first re-construction 
of a gravitational lens using the perturbative approach is described in Alard (2009), this 
work offer a general presentation of the methods and algorithms that will be used in this 
paper. The lens SL2S02176-0513 was already modeled by Tu etal (2009) using conventional 
methods. Tu etal (2009) identified 3 components in the lens: the stellar component, the 
dark matter halo associated with the galaxy, and finally an external influence coming from 
a nearby group of galaxies. The surface density of the bright part is modeled as a Sersic 
profile, while the dark component is represented with an elliptical pseudo isothermal sphere. 
The center of the dark component is supposed to be the same as the center of light. The 
external component is represented with a singular isothermal sphere (SIS). The position of 
the SIS is not a free parameter, it is derived from X ray observations (Geach etal 2007). The 
model proposed by Tu etal (2009) is physically motivated, although it has a large number 
of free parameters and the exploration of the parameter space is always a difficult task. The 
minimization surface in the parameter space is complex, with the usual presence of several 
minima's of similar depth resulting in a degeneracy of the solution. Thus the problem in 
itself is not to find a possible solution but to explore the range of solutions consistent with 
the data. The family of solutions have common properties and these common properties are 
the really interesting quantities. The perturbative approach has the advantage to describe 
the lens using a general singular perturbative expansion. A given perturbative expansion 
corresponds to a family of models, and not to a single model. An important asset in this 
approach is that the general properties of the solution are directly related to the physical 
properties of the images. This direct relation indicates that the perturbative reconstruction 
is not degenerated. A general description of the perturbative methods and of its application 
to the inversion of gravitational lenses is available in Alard (2009). 

2. The perturbative approach in gravitational lensing. 

This section will present a brief summary of the singular perturbative approach in gravi- 
tational lensing, for more details, see Alard (2007). The perturbation is singular since the 
un-perturbed solution is a circle with an infinite number of points, while the perturbed 
solution has only a finite number of points. It is important to note that this perturbative 
approach is not related to conventional regular perturbative method already explored in 
gravitational lensing. For an example of regular perturbative theory in gravitational lens- 
ing, see for instance Vegetti & Koopmans (2009). The un-perturbed situation is represented 
by a circular lens with potential 4>o{r), and a point source with null impact parameter. In 
the perturbed situation the source has an impact parameter rs, and the lens is perturbed 
by the non-circular potential ip(r,9). Both perturbations, rs and tp are assumed to be of 
the same order e: 



[(f) = (j> + e ip 

Similar ideas were explored by Blandford & Kovner (1998), the main difference with Alard 
(2007) is the derivation of analytical equations describing image formation. These equations 
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relate directly the properties of the images to the local potential; the derivation of these 
equations will be presented now. Note that for convenience the unit of the coordinate 
system is the critical radius. Thus by definition in this coordinate system the critical radius 
is situated at r — 1. As a consequence, in the continuation of this work all distances and 
their associated quantities (errors,..,) will be expressed in unit of the critical radius. In 
response to this perturbation the radial position of the image will be shifted by an amount 
dr with respect to the un-perturbed point on the critical circle. The new radial position of 
the image is: 

r = 1 + e dr (2) 

The response dr to the perturbation can be evaluated by solving the lens equation in the 
perturbative regime, leading to the perturbative lens equation: 



/ j r \ dfo 

r s = (« 2 dr - fi) u r - — u e 



(3) 



and k 2 = 1 



d 2 4> 
dr 2 



(r=l) 



This equation corresponds to Eq. (8) in Alard (2007). Considering that the source 
has a mean impact parameter ro, the position in the source plane may be re- written: 
rs = rs + ro- Assuming Cartesian coordinates (xq, yo) for the vector r Eq. breads: 



( K2 dr - A) u r - *k u g 



(4) 



For a circular source with radius Rq, the perturbative response dr takes the simple following 
form (Alard 2007, Eq. 12): 



dr = — 

k 2 



\ 



R2 [ d f° 



(5) 



To conclude it is important to note that Eq. ( [3]) depends on k 2 . However this variable 



can be eliminated from Eq. ( [3|) by re-normalizing the fields: /„ 



_ fn 



, and the source 



plane coordinates, r s = ^ (mass sheet degeneracy). These re- normalized variables will 
be adopted in the continuation of this work. The re-normalization is equivalent to solving 
Eq.'s ( [3]) and ( [5]) for K2 = 1. The variable k 2 will re-appear when the re- normalized 
quantities are converted to the original quantities. 



3. Pre-processing of the HST data. 

The gravitational lens SL2S02176-0513 (Cabanac etal. 2007, SL2S public domain) was 
observed by HST in 3 spectral domain, F475W, F606W, F814W, with an exposure time of 
400 sec. Cosmics cleaning and image re-interpolation to a common grid are performed using 
software's from the ISIS package (Alard 2000). Considering that in the radial direction the 
arc size is only of a few pixels, to facilitate the numerical calculations (convolution with 
the PSf for instance) it is useful to work on a finer grid. The images were re-mapped to 
a grid with a pixel size smaller by a factor of 2, which is small enough to avoid aliasing 
problems. But note that this re-sampling is only a numerical convenience, and that the 
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computing of statistical estimates is performed using the orginal data, and not the re- 
sampled and interpolated (correlated) data. Note also that at this level, no astrometric 
registration is performed, and that the system of coordinate is the system of the initial 
HST image. However the final result will presented in a proper astrometric system (North 
up, and East left). Note also that there is a large number of cosmics in these images, and 
the local density of cosmics is sometime very high. When too many nearby cosmics are 
detected within a small area, the area is flagged, and will be considered cautiously, or even 
rejected. Finally, the 3 cleaned and re-centered images are stacked, and the background is 
subtracted to produce a reference image of the arc system. Two color illustrations of the 
arc system are also provided in Fig.'s (H|)). 




Fig. 1. Color image for SL2S02176-0513 reconstructed from 3 noise filtered HST images in 
3 bands. 

4. Lens inversion. 

The lens SL2S02176-0513 forms a system of image resembling a near cusp configuration, 
with a large arc on the left side of the deflector and a small counter image on the other side 
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(see Fig. [1}. The 4 bright parts in the images have similar brightness and are probably the 
image of a unique area in the source (see Fig. Tu etal (2009) have associated another 
fainter bright detail with these features (see Fig. [J). But we will see that this association 
is unlikely since this detail has a slightly different color and additionally would require a 
complicated field shape to be associated with the other 4 bright features. 




Fig. 2. The red contours indicates the 4 bright areas visible in the image of the source 
formed by the lens. The flux in these areas is similar which suggests that these area corre- 
sponds to the same part of the source plane. The yellow contour indicates a detail associated 
with the 4 other bright details by Tu etal (2009). Note that this detail is fainter and most 
importantly has a slightly different color which makes this association unlikely. 



4.1. Estimation of the critical circle. 

The critical circle is estimated by fitting a circle to the mean position of the images. The 
center and radius of the circle are adjusted. The non-linear adjustment procedure starts 
from a circle centered on the small galaxy at the center of the image. The initial estimate for 
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the circle radius is the mean distance of the images to the center. A non linear optimization 
shows that the best center is close to the initial guess and that the optimal radius is close 
to 30.5 pixels, with a pixel size of 0.049" the critical circle radius Rc is: 1.49". Note that 
the final result will not depend upon the particular choice of a given critical circle. Taking 
another circle close to this one would change a little the estimation of the perturbative 
fields, but the total background plus perturbation would remain the same. 

4.2. General properties of the solution. 

An approximate solution will be derived from the properties of the circular source solution. 
To estimate the circular source solution, the outer contour of the bright area in the source 
will be approximated with a circle. Note that since the bright areas of the image are 
composed of 2 nearby bright spots, the corresponding source region is also composed of 2 
bright features. In the circular approximation, we will consider the circular envelope of these 
2 bright spots. It would have been possible to consider the bright spots one by one, but the 
determination of the local parameters would have been less accurate. This approximation of 
the source may not be very accurate, the typical error on the field estimation will be on the 
order of the source outer contour deviations from circularity. However this approximation 
is sufficient to derive the general properties of the solution, and in particular the behavior 
of the ^ field near the minimum. For circular sources the perturbative fields are directly 
related to the data (see Sec. 3, and Eq. ([5]). The field /i corresponds to the mean radial 
position of the image, while the reconstruction of the field ^ is more complicated and will 
be explained in details in this section. Note that the direct relation between the perturbative 
fields and the data is possible only in bright areas, and that as a consequence ^ and /i 
will have to be interpolated in dark areas. The derivation of the general properties of 
^fy- require some specific guess of the local behavior of this field in the region of image 
formation. Images form in minimum of | | , the local behavior of near the minimum is 
of 2 types: linear behavior for small images , for large images (caustics) ^fe behaves like an 
higher order polynomial (see Alard 2009 for more details). All images of the source bright 
region are small (see Fig [2]), as a consequence a local linear model will be adopted. The 
parameters of the local linear model are estimated using Eq. ( [5]) . Given the values of the 4 
local slopes and the constraint that no image are formed in dark areas ( ^fe > source radius) , 
the general shape of the field is reconstructed (see Fig [3]). Note that the sign of the slope 
for two consecutive images has to be different to avoid a crossing of the zero line between 
the images and the formation of an additional image. Note also that at this level the sign 
of ^jfe is degenerate. Implicitly, since only one feature is seen in the upper right bight area, 
the current model assumes that the 2 bight spots visible in the other images are merged 
here. This is the simplest model consistent with the data. Another possibility is to assume 
that the 2 images are not merged and that the yellow contour in Fig. ( [5]) corresponds to 
the image of the second bright spot, as considered by Tu etal. (2009). However since the 
separation between the images is large the slope of the local model would have to be small. 
To be consistent with this small local slope and the presence of a dark area between the 2 
images, a higher order Fourier expansion of ^ would be required (see Fig. [4|). Additionally 
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the color diagrams of the 4 bright spots (red contours in Fig. ^ are very similar, but the 
color diagram of the fifth bright spot (yellow contour in Fig. [5]) is significantly different 
(see Fig. [SJ- This slight but significant difference in color shows that the fifth bright spot is 
not associated with the same source area as the 4 other bight spots. This is a confirmation 
that the solution proposed in Tu et al (2009) is not consistent with the data, and that that 
the simplest perturbative model (Fig. [3]) predicts the correct association of the images. 




Fig. 3. 




5. Fitting of the perturbative fields. 

A qualitative evaluation of the perturbative fields has been performed in Sec. 14. 21 we will 
turn now to a quantitative evaluation of ^fe . The quantitative evaluation will be conducted 
in two steps, first an approximate guess will be evaluated using the circular source solution, 
and second this guess will be refined by fitting the data. 

5.1. First numerical guess. 

As explained in Section 14.21 the field f± in the circular source approximation is directly 
related to the mean image radial position. To reconstruct this field a Fourier series of 
order 3 is fitted to the mean image position. Fitting a higher order Fourier series does not 
improve significantly the fit, and a lower order expansion is not a good match to the data. 
Let's now turn to the estimation of Sfe. The solution presented in the former section (see 
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Fig. 5. The raw HST flux in each photometric bands normalized by the total flux in the 
3 bands. The red diagram represents the color of the fifth bright area (yellow contour in 
Fig. [2]), the other diagrams represent the color of the 4 other bright areas (red contours 
in Fig. [5]). The color diagram of the fifth bright area is different at the 4cr levels from the 
mean diagram of the 4 other bright areas. 

Fig. [J) is evaluated using the following piecewise numerical model: first order polynomials 
in each area where an image of the source is present, second order polynomials in dark 
areas. Eq [5] shows that the image angular size is given by the condition =fe = R Q . Thus, 
for a linear model, k^A8 = Rq, with A9 the angular size of the image, and k the slope 
of the linear model. The coefficients of the second order polynomials in the dark areas 
are evaluated using continuity conditions, the constraint that no images should be formed 
in dark areas > i?o), and the constraint that the solution should be as smooth as 
possible. The piecewise polynomial solution is fitted with an order 3 Fourier series. The 
Fourier expansion will be used as a first guess for the final fitting of the perturbative fields. 
The numerical guess are presented in Fig O More details about the reconstruction of the 
numerical guess are availbale in Alard (2009), Sec. 4.3 and 4.4. 

5.2. Final fitting of the fields. 

The final fitting of the fields is conducted by reconstructing the images of the source and 
comparing these images with the HST data by chi-square estimation. The reconstruction of 
the images for a given set of perturbative fields will be performed using the Warren & Dye 
(2003) method. The source is represented with a linear combination of basis functions. The 
basis functions are identical to the functionals used for image subtraction (Alard Lupton 
1998, Alard 2000). The images of each basis function is reconstructed and convolved with 
the HST PSF. The HST PSF is estimated using the Tiny Tim software (Krist 1995). 
Finally, the convolved images of the basis function are fitted to the HST data using a 
linear least-square method. The quality of the fit is evaluated by chi-square estimation. 
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Fig. 6. 

This fitting of the data for a given set of perturbative fields is conducted iteratively. The 
inital estimation starts from the numerical guess presented in the former section, leading 
to an initial chi-square estimation. This initial guess is improved iteratively using Nelder 
& Mead (1965) simplex method until convergence is reached. The final result is presented 
in Fig m 

5.3. Noise. 

5.3.1. Chi-square of the fit. 

The model and the data are compared in a small region around the arcs with total number 
of pixels N. The Poisson weighted difference Ri between the model at pixel i, Mi and the 
HST data Di is very close to a Gaussian (see Fig. [3]). Considering that the model has Np 
parameters, the 13 parameters of the perturbative expansion , and the 57 source model 
parameters, the chi-square is, X2/dof — n-n p Si — 1-33. Changing a little the size of 
the area, either reducing it by moving closer to the center of the arcs, or enlarging the area 
does not significantly change the chi-square value. 

5.3.2. Errors on the reconstruction of the perturbative fields. 

Alard (2009) investigated the errors on the reconstruction of the perturbative fields and 
found that the errors due to the Poissonian noise were negligible. The re-construction 
noise is dominated by the errors introduced by the first order approximation of the lens 
equation (see Alard 2009, Sec. 4.6). The results for SL2S02176-0513 are very similar, the 
mean error on the Fourier coefficients of the perturbative expansion due to the Poisson 
noise is: up =~ 0.25 10 -3 (in unit of the critical radius). The error due to the first order 
approximation a a is evaluated by comparing the position of the bright details in the HST 
image and the model reconstruction. The measurement procedure is described in Alard 
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(2009) Sec 4.5.2. The numerical estimation of the error on the position between model and 
data is : ~ 0.6%Rc- This positional error is identical to the error on the perturbative fields 
(see Alard 2009, Sec 4.6.2). Assuming that the errors on the Fourier coefficient are similar 
the error is a a — 2.5e — 03 which is much larger that the error due to the Poisson noise 
op. As a consequence the Poissonian noise will be neglected. Following Alard (2009) Sec 
4.6.3, the error on the reconstruction of the potential iso-contours, S(j)i SO is: 



op 



\ 



1 

^ — ~ 0.0026 

3=1 3 



6. Properties of the lens. 

The perturbative expansion is related directly to the properties of the lens potential. In 
particular Eq. (C.l) in Alard (2009) relates /o to the potential iso-contours. 
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Fig. 8. 

6.1. Potential iso-contours. 

The solution of the perturbative lens equation (Eq [5j) depends on /j = fi+xocos6+yosin6, 
i = 0,1, and Xq, yo, the source impact parameters. The potential iso-contour is related 
to fo by the following equation: : dr = — f (see Alard (2009), Eq C.l). Since only f is 
known, and since the impact parameters are unknown the first order Fourier terms of fo are 
unknown. The first order Fourier terms corresponds to the centering of the potential and do 
not affect the shape of the potential iso-contours (see Alard 2009, Sec 6) .As a consequence 
the shape of the potential iso-contours will be computed using the Fourier expansion of fo 
from order 2 to the maximum order (order 3) . However, the inner potential iso-contour do 
not depend on the impact parameters, as a consequence the relevant first order terms are 
meaningfull. But the corresponding terms are very small, indicating that the center of the 
inner potential iso-contour is very close to the centre of light. An important asset of the 
singular perturbative approach is that the inner and outer contribution to the potential 
can be separated. The potential generated by the lens projected density inside the critical 
circle corresponds to the coefficients a n and b n in the expansion of the potential (see Eq's 
B.l and B.2 in Alard 2009). The potential generated by the the projected density outside 
the critical circle corresponds to c n and d n . Eq B.3 in Alard (2009) relates the Fourier 
expansion of the perturbative fields to a n , b ni c n and d n . As a consequence it is possible 
to reconstruct the potential corresponding to the projected density within and outside the 
critical circle. The potential near the critical circle is dominated by the outer contribution 
(see Fig [9]). Furthermore, the shape of the inner and outer are not the same. Both the 
ellipticity and orientation of the inner and outer contribution are very different. Using the 
second order Fourier terms it is straightforward to estimate the elliptical parameters for the 
inner and outer potential. The elliptical contour equation in a coordinate system aligned 
with the ellipse axis reads: 

(l- V )x 2 + (l + V )y 2 = B 2 o 
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To first order in 77: 




Considering that r = 1 + dr, in a general coordinate system where the mis- alignment 
between the ellipse main axis and the abscissa axis is t/>: 



An identification with the potential iso-contour equation dr = — f leads to, inner contri- 
bution: m = 0.028±0.005 and ipi = -31.7° ±5°, outer contribution: -q = 0.121±0.005 and 
ip a = 9.79° ± 1°. Note that at first order in 77, the ellipticity is equal to 77. It is important 
to remind that as a consequence of the mass-sheet degeneracy the scaling factor K2 is un- 
known (see Eq's [Hand 0. In practice, experiments with dark matter halo's extracted form 
numerical simulations shows that K2 — 1 (Peirani etal 2008). Note that K2 is local quantity 
and that K2 — 1 does not require that the background potential is fully isothermal. In the 
continuation of this paper it will be assumed that K2 = 1, but that in practice K2 may be 
slightly different from unity, which would imply a re-scaling of the perturbative functionals. 
However this re-scaling would affect the ellipticity, but not other properties like the ellipse 
inclination angle. 

6.2. Comparison to observational properties of the lens. 

It is interesting to compare the inner potential iso-contour with the iso-contours of the 
galaxy situated at the center of the lens. To estimate the photometric properties of the 
lens a Sersic profile was fitted to the light distribution. A similar procedure was already 
performed by Tu etal (2009). The surface mass density profile is identical to Tu etal (2009), 
Eq. 1, except that the quantity ~f /n is obviously degenerate and was replaced with a single 
parameter. The iso-contour of this elliptical profile depends on 77 (the ellipticity) and on 
the inclination angle of the ellipse main axis 4>. The fitting of the Sersic profile is conducted 
by minimizing the Poisson weighted residual between the Sersic model and a compilation 
of the 3 HST images. Note that the Sersic model is convolved with the HST PSF, and that 
the non-linear minimization is performed using Nelder and Mead (1965) simplex method. 
The results of this fitting procedure is: r\ = 0.207 ± 0.014 and ip = -33.2° ± 2.4°. Tu etal 
(2009) found: 77 = 0.14 ± 0.06 and ip = 12± 6°. The error bar were evaluated using Monte- 
Carlo simulation (Poisson noise was added to the images). The ellipticity 77 = 0.207 is 
compatible with Tu etal (2009), but the consistency of the ellipse inclination is quite poor. 
The distance between the two measurements is larger than 3 a. Tu etal (2009) suggests 
that the errors may be dominated by systematic effects. To investigate this issue the fit 
was performed separately on the 3 HST images. The results are the following (from blue 
to red photometric band): (-35.1°, -44.9°, -25.1°), with average 35°, and scatter 9.9°. 
Note that estimating the scatter with such a small number of measurements is not very 
accurate. Considering the averaged value of the 3 measurements, no particular systematic 
is detected in this analysis. The reason for the different inclination angle found by Tu etal 
(2009) is unclear. Considering the photometric properties found in the present work, the 
comparison with the iso-contours of the inner potential reads: difference in inclination, 




(6) 
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Fig. 9. The potential iso-contours, the black curve corresponds to the potential generated 
by the whole mass distribution, the blue curve corresponds to the potential generated by the 
outer projected mass distribution and the red curve corresponds the potential generated by 
the projected mass distribution inside the critical circle. Note that the first order Fourier 
terms have not been included in the expansion (due to a degeneracy with the impact 
parameters) . As a consequence the position of the contour center is unknown for the outer 
distribution. In this plot the direction of the vertical axis is parallel to the North (North is 
up and East is left). To improve the readibility of the plot, the amplitude of the non-circular 
distrosion have been exagerated (by a factor of 5) . 

1.5°±5.5°, ellipticities, 0.207±0.014 (galaxy) and 0.028±0.005. Obviously the inclination 
angle are compatible, which suggests that the dark matter halo surrounding the galaxy has 
a similar orientation. The comparison of the ellipticities is model dependent. Assuming the 
following toy model: a dominant isothermal dark halo with ellipticity rjd, to first order in 
r\ the ellipticity of the potential iso-contour is: 3rjd- Thus if the ellipticity of the dark halo 
is the same as the ellipticity of the light distribution, the potential iso-contour ellipticity 
would be: 77 = 0.069 ± 0.005. The difference in ellipticity between the inner potential iso- 
contour is: 0.04 ± 0.007. Thus the ellipticities are not consistent, and this simple model 
suggests that ellipticity of the dark halo is probably smaller smaller that the ellipticity of 
the light distribution. Let's turn now to the contribution of the outer projected distribution 
to the potential. Tu etal (2009) pointed out that the central galaxy is close to a group of 
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galaxies observed by Geach etal (2007). The great difference between the outer and inner 
contribution is due to the influence of this group on the outer component. A simple model 
is to consider two isothermal components for the central galaxy dark matter halo and for 
the mass of the group. Let assume that these 2 isothermal components have respectively 
a velocity dispersion uo and oi, the potential for an isothermal system reads: <j) = Co<x 2 r, 
where Co is a constant. The central galaxy is very close to the center of the coordinate 
system, thus 

4>o = c alr 

in a coordinate system where the center of the group is on the abscissa axis at a distance 
rilll, the potential of the group reads: 

0i = c af\r - n| 



The critical condition implies that = 1, thus Co = -tt, and as a consequence: 

L J r— 1 CT ° 

<t>i = — ||r — r x | 
a o 

To estimate the contribution of the group to the field /o, we will evaluate /q = 4>i(r, 9) at 
r = 1. To second order in — : 



cf>l (r, 9) = -j n + r cos 9 - — 

2 

the coefficient of the second order Fourier term of fn at r = 1 is: — . i . An identification 



with the iso-contour equation (Eq. [SJ, shows that 



CO 

In the former section it was estimated that rj = 0.121, and adopting the position of the 
center of the group proposed by Tu etal (2009), n — 36 (in units of the critical radius). As 
consequence, 21 ==~ 2.95. Geach etal estimate that o\ = 520 ± 120km/s from the a — Tx 
relation and <T\ = 770 ± 170km/s from the X ray temperature. Tu etal. (2009) estimate 
from photometric data that ao = 218^28- Assuming that the mass is dominated by the 
dark halo, respectively: = 2.38 ± 0.72, and 2a. = 3.53 ± 1.04. Both value are compatible 
with 2a. ~ 2.95 inferred from the outer contribution to the potential estimated using the 
perturbative approach. Let's now estimate the compatibility between the inclination angles. 
The group center is situated in a direction at 9q — 16° from the abscissa axis. Since the 
uncertainty on the group center position is 10 arc-sec (Tu etal 2009) this translate in 
an error on Qq of 10.3°. This angle should match the inclination angle of the elliptical 
contour of the outer potential which is: 9.79° ± 1°. Considering the errors, these 2 angles 
are consistent. 



7. Discussion 

The re-construction of the lens SL2S02176-0513 using the perturbative method indicates 
that the gravitational field in this lens is dominated by an outer contribution due to a nearby 
group of galaxies. The separation of the inner and outer potential contribution using the 
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perturbative method is free of any assumptions. Furthermore, the perturbative analysis 
shows that the inner potential contribution is small but statistically significant. The orien- 
tation of the inner potential iso-contour is consistent with the shape of the central galaxy 
luminosity contours, which illustrates the accuracy of the perturbative re- const ruction. 
The outer potential iso-contours reconstructed using the perturbative method are consis- 
tent with the perturbation introduced by the group of galaxies. Both the ellipticity and 
orientation match the group properties. The analysis of this gravitational lens demonstrates 
the ability of the perturbative approach to separate the inner and outer contribution to 
the potential without making any particular assumptions. This lens is very special since 
the outer contribution to the potential is well identified and dominates other contribution 
to the outer potential. In general for other lenses the identification of the perturbators sit- 
uated outside the critical radius is not obvious and using conventional methods the ability 
to re-construct the proper lens potential is compromised. Since the perturbative approach 
does not contain such flaws, and that the inner and outer contribution to the potential 
can be separated it is clear that an un-biased re- const ruction of the lens properties can be 
performed. This capability of the perturbative method to re-construct the inner properties 
of the lens are especially useful to interpret complex gravitational lenses (see Alard 2009). 
Such an analysis will be also very useful in the statistical analysis of the perturbation 
introduced by sub-structures in dark matter halo's (see Alard 2008, and Peirani etal 2008). 
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